function [P] = Crossover(P)
    
    for k = 1:rand(1) * length(P)
        i = ceil(rand(1)*length(P));
        j = ceil(rand(1)*length(P));
        [P(i), P(j)] = Crs_Ind(P(i),P(j));
    end

%     for i = 1:length(P)
%         for j = i+1: length(P)
%             if rand(1) > 0.8
%                 [P(i), P(j)] = Crs_Ind(P(i),P(j));
%             end
%         end
%     end
end

